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Abstract. In this paper we study properties of numbers K l n of connected components 
of bifurcation diagrams for multiboundary singularities B l n . These numbers generalize 
classic Bernoulli-Euler numbers. We prove a recurrent relation on the numbers K n . As 
it was known before, K\ is (n+l)-th Bernoulli-Euler number, this gives us a necessary 
boundary condition to calculate K l n . We also find the generating functions for K l n with 
small fixed I and write partial differential equations for the general case. The recurrent 
relations lead to numerous relations between Bernoulli-Euler numbers. We show them 
in the last section of the paper. 
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1. Introduction 

Like binomial coefficients and Fibonacci numbers Bernoulli-Euler numbers K n are 
widely used in many different branches of mathematics. We define Bernoulli-Euler num- 
bers as Taylor series coefficients for the function sec + tan, namely 



sec t + tant = K n — -. 

' J n I 



n=0 



There exists an equivalent definition of Bernoulli-Euler numbers by means of a "classi- 
cal triangle", one can construct K n by analogy with finding binomial coefficients using 
Pascal triangle. We refer the reader to [1] for the detailed description of the triangle for 
Bernoulli-Euler numbers. For arithmetical properties of Bernoulli-Euler numbers we refer 
to papers [1], [4], and [5]). 



The work is partially supported by RFBR NSh-4719. 2006.1 grant, by RFBR 05-01-02805-CNRSL.a 
grant, NWO-DIAMANT 613.009.001 grant, and by RFBR grant 05-01-01012a. 
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We study one of the geometrical aspects of the numbers K n in singularity theory As it 
was shown by V. I. Arnold in [2], the combinatorics of the components of the set of very 
nice M-morsifications for the boundary singularities of type B n is closely related to the 
combinatorics of the corresponding Springer cones. In [1] V. I. Arnold also proved that 
the numbers of connected components coincides with Bernoulli-Euler numbers. 

In this paper we give proofs of theorems announced by author in [3]. We deal with a 
generalization of boundary singularities B n for the functions on the line to the case of 
multiboundary singularities B l n with boundaries consisting of I points. The numbers of 
connected components of the set of very nice M-morsifications for B l n (denoted by K l n ) 
is in is turn a natural generalization of Bernoulli-Euler numbers. In particular, we prove 
the recurrent relation on the numbers K l n : 

K + -\ = K-nlK l -\ 

Note also that the numbers K l n also enumerate certain strata of singularities of A 2 i+ n -i, 
see Corollary 3.4. 

This work is organized as follows. We start in Section 2 with necessary notions and 
definition, in particular we give definitions of singularities of type B l k . In Section 3 we 
formulate and prove the main theorem on recurrence relation for numbers K l n . Further 
in Section 4 we study the case of small number of boundary points I. We give explicit 
formulae for the numbers K l n in this case. Finally, in Section 5 we show general expressions 
for K l n in terms of Bernoulli-Euler numbers. 

The author is grateful to V. I. Arnold, V. M. Zakalyukin, S. K. Lando, and B. Z. Shapiro 
for useful remarks and help in the realization of this work, and Mathematisch Instituut 
of Universiteit Leiden for the hospitality and excellent working conditions. 

2. Definition of B l n singularities 

Let N m be a smooth manifold, U m be a smooth manifold with smooth boundary, and 
7r : U m — > N m — an immersion, i.e. a smooth embedding with nondegenerate Jacobian 
matrix at any point. We also suppose that the number of preimages of any point of N m is 
bounded from above by some constant. Let / be a smooth real function on the manifold 
N m . Denote by / the lifting / o ir : U m -> K. 

Let the preimage of the point x in N m consists of I boundary points of U m and any 
number of nonboundary points of U m . Let also the collection of / images of tangent 
hyperplanes to the boundary points projecting to collection of hyperplanes in 

general position: the intersection of any s (for s < I) hyperplanes is a (m— s)-dimensional 
plane. 

Definition 2.1. We say that the function / has a multiboundary singularity of type B l n 
at the preimage 7r _1 (a;o), if / has a singularity A n _i at point x . In addition it is required 
that the kernel of the Hessian matrix at xq is transversal to the projection of any tangent 
hyperplane at the boundary point in the preimage 7r _1 (:ro). 

Let us remind a definition of a very nice M-morsification of the boundary singularity 
B^ of the paper [2]. 
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Consider the space R M_1 of real polynomials with zero constant turm 

x» + Aix^ 1 + . . . + A M _xx 
on the real line with fixed "boundary" x — 0. 

Definition 2.2. A very nice M-morsification of the boundary singularity is a polyno- 
mial of this family whose all /i— 1 critical points are real and the values at these points 
are distinct and nonzero (i.e. also distinct to the value at the boundary x — 0). 

Let us generalize Definition 2.2 to the case of multiboundary singularities of type B l n . 
Consider the Cartesian product of the spaces of real polynomials (M" - " 1 ) 

x n + \ 2 x n ~ 2 + ■■■ + \ n _ lX 

and the space M. 1 of the space of "boundary values" at the points x = bi, for % — 1, . . . , I. 
We call the points bi the boundary points. 

Definition 2.3. A very nice M-morsification of the multiboundary singularity B l n is a 
degree n polynomial with / marked boundary points, whose all n—\ critical points are 
real, and the values at critical points and at I boundary values are all pairwise distinct. 

Notice that we enumerate all boundary points, otherwise we should consider the factor 
of the space ~R l by the group of all permutations of the coordinates. We enumerate the 
boundary points since they correspond to distinct branches, whose permutations do not 
make sense here. 

Definition 2.4. The M-domain is the closed subset of the space of all polynomials x n + 
\ 2 x n ~ 2 + • • • + \ n ~i x , consisting of polynomials which critical points are real. 

The set of very nice M-morsifications is an open subset of R™ -1 x M 1 , and its closure is 
(M-domain) x M!. The set of very nice M-morsifications splits in connected components 
by the bifurcation diagram consisting of five hypersurfaces. Three of these hypersurfaces 
come from the case of boundary singularities of type B n (see also in [2]): 

(a) : the boundary caustic consisting of functions with a boundary critical point; 

(b) : the ordinary Maxwell stratum consisting of functions with equal critical values 
at different points; 

(c) : the boundary Maxwell stratum consisting of functions having some value at the 
boundary being equal to some critical value (the corresponding critical point is 
not at the boundary). 

Notice that in our notation B n = B\, moreover the boundary point does not fixed at zero. 
So the definition of the very nice M-morsification of the boundary singularity B n in the 
paper [2] is a special case of Definition 2.3. 

In the case of multiboundary immersion singularities of type B l n for / > 2 we have 
additional two hypersurfaces: 

(d) : the double boundary caustic consisting of functions with some double boundary 
point; 

(e) : the double boundary Maxwell stratum consisting of functions with equal values 
at different boundary points. 
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3. Recurrent relation on numbers Ki 



We use K n to denote Bernoulli-Euler numbers (let us list some first numbers K n starting 
from n = 0: 1,1,1,2, 5, 16, 61, . . .). Denote by K l n the number of connected components of 
the set of very nice M-morsifications of the singularity B l n . The Bernoulli-Euler numbers 
are boundary conditions for the numbers K l n , namely K® = K n _i and K\ = K n+1 (see 



We start the proof with two lemmas. In these lemmas we use the following notation. 

Denote by L l n the number of connected components of the boundary caustic in the 
complement to the union of the strata and caustics of codimension 2. Here we suppose 
that the function has one critical point at the projection of a boundary point, all the rest 
critical and boundary points are distinct and the values at them are also distinct. 

Lemma 3.2. Let n > 2 and I > 1. Then the following identity holds. 



Proof. We take any very nice morsification with /— 1 boundary points and n—1 critical 
points. Consider the action of the permutation group on the boundary points. This action 
naturally defines (I— 1)! very nice M-morsifications in different connected components but 
with the same set of boundary points and the same polynomial. We put a new Z-th 
boundary point to one of n—1 critical points. The action of the permutation group on the 
boundary points now defines /! distinct M-morsifications with the same set of boundary 
points and the same polynomial. Therefore, we uniquely associate to the collection of 
(I — 1)! old connected components the collection of (/ — 1)! new connected components. 
This implies the statement of Lemma 3.2. □ 

Further we find the number of connected components of very nice M-morsifications for 
which one of the boundary values is either greater than all the critical values or less than 
all the critical values, we denote this number by K l n . Note that we count twice connected 
components of all M-morsifications for which all the critical values are contained in the 
interval with endpoints at two boundary values. 

Lemma 3.3. Let n > 2 and I > 1. Then the following identity holds. 




L l n = l(n - l)K\ 



n 



k; 



2ik: 



i-i 



n 



n 
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Proof. Let n be even and I > 1. Then any M-morsification has two branches tending to 
plus infinity. Consider the maximal boundary value. If it is greater than all the critical 
values than it can be attained only at a point corresponding to two branches described 
above. Hence, for any very nice M-morsification with /— 1 boundary point we have exactly 
I distinct M-morsifications for which one of the boundary points bi where 1 < i < I is on 
the right branch, and the value at it is the maximal among all the boundary and critical 
values (here we keep the order of the rest points). The same reasoning is valid for the 
boundary point with maximal value at the left branch. Hence, we obtain K\ n = 21K 1 ^ 
for even n. 

If n is odd, then any M-morsification has two branches, one of them tends to plus infinity, 
and the other tends to minus infinity. By the same reason the statement of Lemma 3.3 
hold for odd n. We remind that we count connected component corresponding to the odd 
case twice. □ 

Proof of Theorem 3.1. Consider a connected component of the set of very nice M- 
morsifications with /+1 boundary points and n— 3 critical points. Take any M-morsifica- 
tion of this component. Consider the last boundary point b l+ i and add to it a 8- 
shaped function concentrated in a small neighborhood of this point such that the new 
M-morsification has a critical point at b l+ i with old boundary value. The value of the 
second critical point that occurs while adding the 5-shaped function is the maximal among 
all critical and boundary values. New M-morsification has now I boundary and n— 1 crit- 
ical points. In the same way we can subtract 5-shaped function. This provides us the 
decomposition of the set of M-morsifications with I boundary and n— 1 critical points 
into couples (corresponding to addition and subtraction of 5-shaped function). For any 
connected component of the set of very nice M-morsifications with l+l boundary and 
n— 3 critical points by the above procedure we bijectively associate a couple of connected 
components of very nice M-morsifications with / boundary and n— 1 critical points. 

Consider an arbitrary connected component of the set of M-morsifications of the bound- 
ary caustic with n— 1 critical point and I boundary points. So, one of the boundary points 
is critical, and in other boundary and critical points the values are pairwise distinct. Con- 
sider the critical boundary point. Let us move the boundary point a little to the right or 
to the left. If in the critical point we have local maximum than we add 5-shaped function 
in this critical point, otherwise we subtract 5-shaped function. The critical value becomes 
then maximal (or minimal respectively). For the boundary point there are two positions 
in correspondence with the direction we have moved it. This gives us to two very nice 
M-morsifications with n— 1 critical points and I boundary points. Hence, we associate 
to each of the components of the boundary caustics a couple of components of very nice 
M-morsifications with maximal or minimal critical values. 

Finally, consider a component of very nice M-morsifications with n— 1 critical points 
and / boundary points with one critical value greater than all boundary values. Consider 
an M-morsification of this component. Let Xi be the point with maximal critical value. 
Let us take critical or boundary points that are the closest to Xi from the left and from 
the right respectively. Choose one of this two points with maximal value. Let us make an 
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operation inverse to the adding 5-shaped function that "pulls down" the maximal critical 
value to the level of the neighbor chosen point. The result can be of two types. The first 
possibility is if we obtain a function with double critical point. In this case we replace 
a double point by a new boundary point b\ +1 . The second possibility is if we obtain a 
function having one critical point coinciding with one boundary point. 

We act in the same way for the case of a critical value being less then a;; boundary 
values. 

The above observations prove the following identity: 

2K l n -K l n = 2L l n + 2K l +\ 
Now we apply Lemmas 3.2 and 3.3: 

2K l n - 21K 1 ' 1 = 2l(n - \)K l ~ x + 2K l n + _\. 

Therefore, 

K l n + _\ = K l n -nlK l ~\ 

This concludes the proof of Theorem 3.1. □ 

In conclusion of this section we indicate the relation between K l n and the numbers 
of connected components of special strata of bifurcation diagram of critical points and 
critical values for degree 2l+n polynomials (i.e. of A 2 i +n ~i)- 

Corollary 3.4. Consider open strata of caustic of the singularity A 2 i+n-i corresponding to 
polynomials of M-domain with l-couples of double critical points with all distinct critical 
values at distinct critical points. The number of connected components of such strata 
tsK l Jl\. 

4. Corollaries of Theorem 3.1. Cases of small numbers of boundary 

points 

Let us give explicit formulae for K l n for / < 5. In [1] V. I. Arnold proved that the number 
of connected components of very nice M-morsifications for singularity A n _i equals to the 
n-th Bernoulli-Euler number, namely K(A n _i) = K n _ x - It [1] there is also a formula for 
numbers K(B n ) for multiboundary singularities of type B n : K(B n ) = K n+1 . So, we have: 
= K n - U and K\ = K n+1 . 

Direct calculations lead to the results of the following corollary of Theorem 3.1. 

Corollary 4.1. The following expressions of K l n in terms of Bernoulli-Euler numbers are 
correct. 

K 2 n = K n+3 -(n + 2)K n+l] 
Kl = K n+5 - (3n + 8)K n+3 ; 

K A n = K n+7 - (6n + 20)^ n+5 + 3(n + 2)(n + A)K n+3 ; 

Kl = K n+9 - (lOra + 40)fT n+7 + (15n 2 + UOn + l84)K n+5 
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The exponential generating function for Bernoulli-Euler numbers is the function 

K{t) = tan(t) +sec(f). 

Let us show exponential generating functions for the numbers B l n for / < 4. Note that in 
the cases of I = and I = 1 one can propose the exponential generating functions to be 
K(t), still in our notation we have 

K (t) = J K(t)dt = - ln(cos(t)) + ln(tan(| + f )) + C and 

KAt) = K'it) = ^f- = ±±^# = "j — ^777 
lv - J W 2 cos^ t 1— sm(t) 

respectively. 

Corollary 4.2. The exponential generating functions for the cases I = 2,3,4 are as 
follows: 

K 2 (t) = K"'(t) - (tK(t)y = 3s ffr;g;f) ; 

K 3 (t) = (K" - 3tK' + K)"' = (1 _ s . 3 n(t))3 (sin(t)(3sin(t) + 7) - 3tcos(t)(5 + sin(t))); 
K A {t) = {K"> - UK" + (3t 2 +4)K'- 3tK)^ = 



l-sin(t) (l-sin(t)) 2 V u ^^vjj i ( 1 _ sin ( t )) 

Let 

K(x,y) = J2lTi xl y n 
l\n\ 

be an exponential generation function in two variables. 

Corollary 4.3. The function K(x,y) satisfies the following differential equation 

K X = (1 1x)Kyy Xt/Kyyy. 

B. Z. Shapiro proposed to consider two exponential generating functions in two variables 
R(x, y) and S(x, y) separately for R l n = K l 2n and for S l n = K l 2n _ v This decreases the order 
of the differential equation. 

Corollary 4.4. The functions R(x,y) and S(x,y) satisfy the following differential equa- 
tions: 

R x = (1 - 2x)Ry - 2xyR yy , 
S x = {l-2x)S y -x{2y-l)Syy. 

Finally, let us describe geometrical structure of B\ with schematic "pictures" of fibers 
in the bundle 

7r : (M-domain) x R 2 — > M-domain. 

For any polynomial f(x) we associate a skew-symmetric polynomial f{b\) — /(&2) i n t w0 
variables b\ and b 2 . This polynomial defines the double boundary Maxwell stratum and 
boundary caustics in the fiber. Boundary caustics and and boundary Maxwell stratum are 
defined by vertical and horizontal straight lines and form a rectangular net, in which the 
curve f(bi) = /(&2) is inscribed. It is easy to draw combinatoric pictures of such curves. 
For example, we draw the pictures for the singularities B\ and B\ (see Figures 1 and 2). 
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Notice that the curve fipi) = f(b 2 ) is a union of a straight line 61 = b 2 with some curve 
of degree n—1. For instance, in the case of B\ it is the union of a straight line b\ = b 2 
and an ellipse (see Figure 1). Here a natural problem arises. Describe all combinatoric 
types of such pictures for general n. At present moment we do not know the answer to 
this problem. 

5. Corollaries of Theorem 3.1. Connection with Bernoulli-Euler 

NUMBERS 

In this part we introduce some expressions for the numbers K l n similar to the expressions 
of Corollary 4.1. Further we calculate K l n for negative n satisfying n < —I. These numbers 
lead to nice relations on Bernoulli-Euler numbers. 

Corollary 5.1. The following expression for K l n is true: 

K l n = K n+ 2l-l — 

(, — 2 — 3 ) J ^n+2l-3 + 

1(1 -m~ m - 3)(|n 2 + 2ttln + {l+l) s~ 2) )K n+2l _ b + 

E ( L^nl-iy. {i-2k+2)\ E iPk,d{l)n d j) K n+2(l _ k)+l ) , 

k=4 y y d=0 ' ' 

where Pk,d{ x ) is a polynomial of degree k—d—1 with constant coefficients depending on k 
and d. 

The idea of the proof is based on the induction on /. 

Still we do not know the explicit formulae for the polynomials Pk,d(x), nevertheless there 
exists a recursive method to find the coefficients of such polynomials. Let us show the 
polynomials for d = k—1, d = k—2, and d = k—3. 

Corollary 5.2. The following hold: 

p k ,k-i(x) = 1; 

Pk,k-2{x) = ^(2x + 4-k); 

p k , k - 3 (x) = (fc -^ fc ' 2) (20x 2 + (72 - 20k)x + (5k 2 - 39k + 64)) . 

Applying formulae of Corollary 5.1 one can obtain many relations for Bernoulli-Euler 
numbers. Let us consider some examples of such relations. Substitute n — 1 and n = 2 
in the formulae of Corollary 5.1, we get the equalities of the following theorem. 

Theorem 5.3. The following relations on Bernoulli-Euler numbers hold: 

k n fc-l X 

k 21 - 1\ = k 2 i - k[ = ( 2 *-i(fc-i)i a-2k+2)\ E PkA l )) K 2(i~k+iy, 

k=2 V d=0 7 

K 2l+1 - 2H\ = K 2l+1 -K\= E ( aJr^m ( ^ +2)! E(p M (02 d ))^- fc+ i )+ i. 

k — 2 d — 
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Figure 1. A fiber of general position of the M-domain for the singularity B\. 




Figure 2. A fiber of general position of the M-domain for the singularity B\. 



Let us use a expressions of Corollary 4.1 to obtain the numbers K l n , for n < 0, and 
n < -I. 
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The numbers K can be interpreted as numbers of connected components of the space 
of degree polynomials with I boundary points with distinct boundary values. Here we 
have only one non-zero entry corresponding to I = 1. The author does not know any 
reasonable interpretation for the numbers for negative n. 

Notice that there are many zeroes in the table. Let us formulate this statement more 
precisely. 

Proposition 5.4. Let n < -I, then K l n = for I > -n and K~ n = (-n-l)L The 
numbers K l = for I > 1 . 

Proof. The proof is based on Theorem 3.1. Let us use the induction on —n. 

K l = K l 2 ~ l - 2(1 - l)K l f 2 = {21 - 2)!! - (21 - 2)(2l - 4)!! = for I - 2 > 0. 
K l _ x = K 1 ^ 1 - (I - l)K[- 2 = (I - 1)! - (/ - 1)(Z - 2)! = for I - 2 > 0. 
K l _ 2 = K 1 ^ 1 = for I - 1 > 2. 

In general case we have 

K = K l - + \ -( n + 2)(l- l)K l ~ + 2 2 for n< -2, 1 - 2 > -2 - n. 

K~ n = K^ 1 - (n + 2)(-n - 1)K^ 2 2 = + ( n + 2)(n + l)(-n - 3)! = (-n - 1)!. 
This proves the statement of the proposition. □ 

We conclude the paper with relations for the Bernoulli-Euler numbers that follows from 
the results of Corollary 5.1 and Proposition 5.4. 

Corollary 5.5. Let n < 0, I > max(l, —n), then 

= K = K M + £ {( J^ k _\ y {l _^ k + 2)] E(PM(0- d ))^ +2 ,- fe)+ i)- 

k=2 K ' V ; d=0 

If I = —n, then 

(/ - 1)! = at = AV, + E ((^(T^EWOI-O'))^). 

fc=2 V J K ' d=0 

In particular for I > 1 we have: 

(-l) fe_1 /' \ \ 

= /^ = K n+2l _ ± + g (( 2fc _ 1(ifc _ 1} , (z _ 2fc + 2) ,P fc ,o(0 j^2(i-* )+ ij ■ 
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